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We use series expansion methods to study ground- and excited-state properties in the helically 
ordered phase of spin-1/2 frustrated antiferromagnets on an anisotropic triangular lattice. We 
calculate the ground state energy, ordering wavevector, sublattice magnetization and one-magnon 
excitation spectrum for parameters relevant to CS2CUCI4 and Cs2CuBr4. Both materials are mod- 
eled in terms of a Heisenberg model with spatially anisotropic exchange constants; for CS2CUCI4 
we also take into account the additional Dzyaloshinskii-Moriya (DM) interaction. We compare our 
results for CS2CUCI4 with unpolarized neutron scattering experiments and find good agreement. 
In particular, the large quantum renormalizations of the one-magnon dispersion are well accounted 
for in our analysis, and inclusion of the DM interaction brings the theoretical predictions for the 
ordering wavevector and the magnon dispersion closer to the experimental results. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

The study of frustrated quantum antiferromagnets is 
central to modern condensed matter physics. Much of the 
interest in these many-body systems stems from the pos- 
sibility, first envisaged by Anderson^ that their quantum 
fluctuations may be so strong as to lead to exotic "spin- 
liquid" ground states characterized by the absence of bro- 
ken symmetries of any kind. Another hallmark signature 
of spin liquids is the existence of deconfined fractionalized 
spin-1/2 excitations, usually referred to as spinons. Un- 
fortunately it has been very hard to find experimental re- 
alizations of such states. Important progress has however 
been made recently with the identification of some ma- 
terials which may have, or at least be close to having, a 
spin-liquid ground state. Two such promising candidates 
are the organic compound k-(BEDT-TTF) 2 Cu 2 (CN) 3 2 
and the layered S — 1/2 frustrated antiferromagnet 
Cs 2 CuCl 4 . 

The magnetic properties of CS2CUCI4 have been ex- 
tensively studied using neutron scattering^ At low tem- 
peratures T < Tjy = 0.62 K long-range helical mag- 
netic order is observed and the low-energy excitation 
spectra contain relatively sharp modes characteristic of 
magnons, the expected Goldstone modes. The magnon 
dispersion does however show very strong renormaliza- 
tions compared to linear spin-wave theory. Furthermore, 
the dominant feature of the spectra is in fact a broad 
continuum, occurring at medium to high energies, which 
carries most of the spectral weight and which persists 
also for temperatures above TV. The interpretation of 
this continuum has recently been the subject of much 
debate^&l&^'I^ 2 .^ 

The most interesting hypothesis, already suggested in 
Ref. [H, is that the continuum is due to two-spinon scat- 
tering resulting from CS2CUCI4 being close to a quan- 



tum phase transition to a two-dimensional spin-liquid 
state. A number of theoretical proposals have been made 
regarding the nature of the spin liquid that might be 
involved in such an unconventional scenario4'£' 7 -' 8 - : 11 i 12 
An alternative hypothesis, that the effects of magnon- 
magnon scattering included within a standard (nonlin- 
ear) spin- wave approach might be able to explain the ex- 
perimental results, was recently explored in Refs. [9| and 
[10 . While qualitative features of the calculated spec- 
tra were found to be similar to the experimental ones, 
the quantitative agreement was however not satisfac- 
tory. In particular, the obtained quantum renormaliza- 
tion of the magnon dispersion (calculated to lowest order 
in magnon-magnon interactions) was significantly lower 
than found experimentally. 

The magnetic interactions in CS2CUCI4 are sufficiently 
weak that the external magnetic field required to fully po- 
larize the spins is experimentally accessible. By measur- 
ing the magnon dispersion relation in this fully-polarized 
state (in which quantum fluctuations are completely sup- 
pressed by the applied field) , it was found 14 - that the lay- 
ers in CS2CUCI4 are nearly decoupled from each other and 
well described by a spin-1/2 triangular-lattice Heisen- 
berg antiferromagnet with exchange constants J and J' 
(defined in Fig. [1]) with J'/ J w 1/3, weakly perturbed 
by an additional intra-layer interaction of Dzyaloshinskii- 
Moriya (DM) type^ of strength D/J w 0.05. Although 
weak, the DM interaction has several important conse- 
quences, one being that the ordering plane of the spins 
coincides with the plane of the triangular lattice because 
the DM interaction makes the latter an "easy" plane. 

Theoretical studies of CS2CUCI4 and models rele- 
vant to it have been carried out using many different 
methods. For the zero-field case, which is our focus 
here, these methods include spin-wave theor y 9 ! 10 : 16 series 
expansions, 17 ^ variational approaches ) 12 ' 18 and various 
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field-theoretical approaches^^ ^ 11 ' 13 ! 19 While most of 
these studies only considered the Heisenberg part of the 
Hamiltonian, some of them also investigated the effects 
of the DM interaction. Rcfs. HI andflol using the nonlin- 
ear spin-wave approach, found that several quantities, in- 
cluding the ordering wavevector and the quantum renor- 
malization of the magnon dispersion, were quite sensitive 
to the DM interaction, and that the reduced spin rota- 
tion symmetry of the Hamiltonian for D ^ also leads to 
a gap in the magnon dispersion at the ordering wavevec- 
tor. Ref. [13 also studied the sublattice magnetization 
as a function of D (for values of J, J' appropriate for 
CS2CUCI4) and found that it vanishes when D/J< 0.008, 
and hence that the DM interaction is crucial to stabilize 
the magnetic order in CS2CUCI4. The same qualitative 
conclusion was reached in Ref. 13 which used a quasi- 
one-dimensional approach, treating the interchain inter- 
actions (J' and D) using a perturbative renormalization 
group analysis; these authors found that in the absence 
of the DM interaction the ground state was dimerized. 
Finally, the easy-plane nature of the Hamiltonian when 
D 7^ plays an essential role in the algebraic vortex- 
liquid scenario proposed for CS2CUCI4 in Ref. [ill : also in 
this study the magnetic order was found to be driven by 
the DM interaction. 

Another interesting material is Cs2CuBr4, which can 
be described in terms of a spin- 1/2 triangular-lattice 
Heisenberg antiferromagnet with J'/ J w 1/2^ Com- 
pared to CS2CUCI4 this material is therefore more frus- 
trated (i.e., closer to the isotropic limit J = J') and fur- 
ther away from the one-dimensional limit J J' . The 
thermodynamic properties of both materials, particularly 
the uniform magnetic susceptibility, have been exten- 
sively studied. There are also a large number of organic 
materials from superconducting families that have insu- 
lating magnetic phases that can be described by a spa- 
tially anisotropic triangular-lattice models Recent high 
temperature series expansion studies of these systems^ 
found that a large class of the organic materials appeared 
close to the isotropic triangular-lattice limit, while the 
inorganic materials were closer to the one-dimensional 
limit. 

In this paper we study spin-1/2 triangular-lattice an- 
tiferromagnets using zero-temperature high-order series 
expansions. We focus on the case J > J' for which, 
according to our previous series expansion study in Ref. 
[l7l the ground state has noncollinear (helical) long-range 
magnetic order and the ordering wave vector varies con- 
tinuously with the model parameters. We calculate 
the ground state energy, ordering wavevector, sublat- 
tice magnetization, and magnon dispersion for values of 
J'/ J relevant to Cs 2 CuCi4 and Cs 2 CuBr 4 . For Cs 2 CuCl4 
we also consider the effects of the DM interaction. We 
compare our results with those of other theoretical ap- 
proaches and with neutron scattering experiments. A few 
of the results discussed in this paper have already been 
briefly presented in Ref. I22I 

The outline of the paper is as follows. In Sec. HT1 we 



describe the model Hamiltonian. In Sec. IIIII we dis- 
cuss the series expansion method used for studying zero- 
temperature properties in the helical phase. In Sees. IIVI 
and [V] we discuss our results for ground state properties 
and the magnon dispersion, respectively. Finally, our 
conclusions are presented in Sec. IVII 
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FIG. 1: Exchange constants J, J' and lattice basis vectors 
Si, 82 on a triangular lattice. The standard notation for axes 
in CS2CUCI4 is also indicated. 



II. MODEL HAMILTONIAN 

The antiferromagnetic spin-1/2 Heisenberg model with 
spatially anisotropic exchange constants on a triangular 
lattice forms the basis for a theoretical description of the 
magnetic properties of both CS2CUCI4 and Cs2CuBr4. 
The Heisenberg Hamiltonian reads 

#h = ^ R ' Sr+S!+S2 + J' Sr ' {Sr+S! + Sr+s 2 )}j 

R 

where Sr is the spin operator on site R, J and J' are 
the antiferromagnetic exchange constants, and S± and 82 
are nearest-neighbor lattice vectors, as defined in Fig. [1] 
When the ground state of has magnetic long-range 
order, the spin expectation value can be written 

(S R ) = M[ni cos(Q • R) + n 2 sin(Q • R)} (2) 

where M is the sublattice magnetization, n\ and n-2 are 
two arbitrary orthogonal unit vectors, and Q = (0, Q, 0) 
is the ordering wavevector with cos(Q6/2) = cost?, where 
q (2q) is the angle between nearest-neighbor spins cou- 
pled by J' (J) and b is the lattice constant in the b direc- 
tion (see Fig. [I]). For classical spins (corresponding to the 
limit of spin S — > 00), M — S and q — arccos(— J'/2 J) 
(for J > J'/2). Eq. ([2]) then implies that the magnetic 
order is helical, with the ordering plane spanned by the 
vectors n\ and ri.2. For quantum spins, quantum fluctu- 
ations could in principle kill the magnetic order. How- 
ever, in one of our earlier series expansion studies^ of 
the model (JTJ) for S — 1/2, the magnetic order appeared 
to be robust as long as J / J' is not too large. The quan- 
tum fluctuations do however renormalize Q (and hence q) 
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away from the classical value and reduce the strength of 
the ordering so that the sublattice magnetization M < S. 

Cs2CuBr4 can be described by the Hamiltonian ([1]) 
with J I J' k. 2. 20 In contrast, Cs 2 CuCl 4 has J = 0.374(5) 
meV and J' = 0.128(5) meV, giving a more anisotropic 
ratio J I J' « 2.92. 14 Furthermore, the spin Hamiltonian 
of CS2CUCI4 also contains a Dzyaloshinskii-Moriya (DM) 
interaction 1 ^ -£?dm of the formii 

Hum = ' Sr x {Sr+Sx + Sr + s 2 ), (3) 

R 

where D lies along the a-direction, perpendicular to the 
plane of the triangular lattice, with D — 0.020(2) meV, 
giving D/J' sa 0.16. Thus the DM interaction is numer- 
ically a relatively small perturbation on the dominant 
Heisenberg Hamiltonian Hjj- Nevertheless, the DM in- 
teraction has several notable consequences. It breaks the 
full SU(2) spin rotation symmetry of the Heisenberg part 
down to U(l) by making the plane of the triangular lat- 
tice an "easy plane" which therefore becomes the order- 
ing plane of the spins. Furthermore, as the DM interac- 
tion can be seen to give rise to a linear coupling to the 
unit vector n\ x n 2 (which points perpendicular to the 
ordering plane), it selects a unique direction for this vec- 
tor, corresponding to a specific chirality or handedness 
of the spin order. In CS2CUCI4 the direction of D in fact 
alternates from layer to layer— and hence the chirality 
alternates as well. 

When discussing CS2CUCI4 in this paper we will con- 
sider both the Hamiltonian Hh + ^dm with D/J' = 0.16 
as well as the more simplified model defined by neglect- 
ing the DM interaction, i.e., D = 0. In CS2CUCI4 there 
is also an antiferromagnetic interlayer exchange interac- 
tion J",— but because this is very small (J" = 0.017(2) 
meV w -Z/45) and because its effectiveness in coupling 
the layers is further reduced by the alternating chirality, 
we choose not to include this interlayer interaction in our 
analysis, thus focussing on a single layer. 

III. SERIES EXPANSIONS IN THE HELICALLY 
ORDERED PHASE 

In this section we discuss some aspects of the series 
expansion analysis of the helically ordered phase. 

Following Ref. [23|, we assume that the spins order in 
the x z plane (which would be the be plane in CS2CUCI4), 
with D pointing in the perpendicular direction. We ro- 
tate all the spins so as to have a ferromagnetic ground 
state, with the resulting Hamiltonian H: 

H = H 1 + JH 2 + H 3 (4) 

where 

Hi = J cos (2q) S?S z n +[J' cos (q)-D sin(q)] ]T S*S| , 

(in) (ij) 

(5) 



H 2 = SfSy + cos (2q)SfSZ + sin (2q)(SfS^ - SfS*), 

(in) 

(6) 

(ij) 

+ [J'sm{q) + Dcos{q)]{S*S* - SfS?). (7) 

where the sum (in) is over nearest-neighbor sites con- 
nected by "horizontal" bonds in Fig. [1] with exchange 
interaction J, and the sum {ij) is over nearest-neighbor 
sites connected by "diagonal" bonds with exchange in- 
teraction J'. 

Next, we introduce the Hamiltonian 

H(X) =H + XV (8) 

where 

i/o = Hi-tE(^- 1 /2), (9) 

i 

V = JH 2 + H 3 +tJ2(S? ~ 1/2). (10) 

i 

The last term of strength t in both Hq and V is a lo- 
cal field term, which can be included to improve conver- 
gence. H(\ = 0) is a ferromagnetic Ising model with two 
degenerate ground states, while H(\ = 1) is the model 
whose properties we are interested in. We use linked- 
cluster methods to develop series expansion in powers of 
A for ground state properties and the magnon excitation 
spectra. 

For J J', the lattice has C 2v symmetry (the symme- 
try operations are the identity, inversion, and reflections 
about the b and c axes) , and the series for the spin-triplet 
excitation spectra has the following form: 

Uj{k x ,k y )/J' = ^ A r ^ C r , m ,n C0S (tT M COS (—^- k v) 
r— m,n 

(11) 

where c Ttm ^ n are series coefficients, m and n are integers, 
representing a hopping over distance (m/2, nV3/2), and 
m + n is a even number. For this case, series for the spin- 
triplet dispersion has been computed to order A 8 , and the 
calculations involve a list of 25022 linked-clusters, up to 
9 sites. The series coefficients c r:in . n for J = 2.92, J' = 1, 
q = 1.64, D = 0.16, and t = 4 are given in TableQH 

For more details we refer to Ref. [23j where series expan- 
sions for the ground state properties and magnon dispe- 
rion of the spin- 1/2 Heisenberg model on the anisotropic 
triangular lattice are discussed at length^ 

IV. GROUND STATE PROPERTIES 

A. Ordering wavevector 

In Fig. [2]we present results for the ground state energy 
per site, as a function of the ordering angle q, for various 
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FIG. 2: The ground state energy per site, as a function of 
the angle q between nearest-neighbor spins along J' bonds, 
obtained from the series expansion around the Ising limit, 
(a) J I J' = 2, relevant for Cs2CuBr4. Here the left and right 
vertical dashed lines denote the value of q obtained from dimer 
expansions [Eq. (16) in Ref. [TtJ, and the classical ground 
state, respectively, (b) J/ J' = 2.92, the ratio appropriate 
for Cs 2 CuCl 4 . The DM interaction (D = 0.16J') lowers the 
ground state energy and shifts the ordering angle from q w 
1.59 to q » 1.64. 



model parameters. The actual, realized, value of q is that 
which minimizes the ground state energy. 

CsiGuBr^ [Fig. HJa)]. The value of q predicted by the 
Ising expansion is seen to be considerably smaller than 
the classical value (right dashed line) but is quite close to 
the value predicted from a dimer expansion (left dashed 
line; Eq. (16) in Ref. G3). 

CsiCuCl^ [Fig. [2^b)]. The ordering wavevector Q 
(equivalently, the ordering angle q) is conventionally ex- 
pressed in terms of the incommensuration e, defined as 
the deviation between Q and the antiferromagnetic (with 
respect to the b direction) wavevector n/b, measured in 
units of 2n/b: e = (Q - n/b)/(2n/b) = q/n - 1/2. For 
D = the series calculation gives q m 1.59 which implies 
e w 0.006, while for D/J' = 0.16, q « 1.64 which gives 



e w 0.022. Thus inclusion of the DM interaction increases 
the incommensuration and brings it much closer to the 
experimental value e = 0.030(2). 14 

In Table U we compare our results for the incommen- 
suration e for the helically ordered phase in CS2CUCI4 
with experiments and with predictions from some other 
theoretical approaches. 



B. Sublattice magnetization 

For Cs 2 CuCl 4 (J/ J' = 2.92 and D/J' = 0.16) we find 
the sublattice magnetization to be M = 0.213(10). In 
this result the main source of error is the uncertainty in 
the angle q. It is interesting to compare this with predic- 
tions from other theoretical approaches (see Table Q| . In 
particular, by using the 1/5 expansion and taking into 
account quantum corrections up to and including order 
1/S 2 to the classical result, Ref. llCj found the slighly 
smaller value M sw 0.18 for the same parameters, and 
also predicted that M = for D = 0, i.e. that quantum 
fluctuations would be strong enough to completely melt 
the magnetically ordered state in the absence of the DM 
interaction. In contrast, our analysis indicates a small 
but nonzero value M ss 0.10 for the magnetization when 
D = 0, although it should be noted that the error bars 
are significant (see Fig. 6 in Ref. 17). We conclude 
that in both approaches the DM interaction is found to 
cause a considerable suppression of quantum fluctuations 
and strengthen the magnetic ordering tendencies. Includ- 
ing the interlayer coupling J" is not expected to change 
the magnitude of the ordered moment by much, as J" is 
rather small (J" / J ~ 45) and also because the chirality 
of the spin ordering is opposite in neighboring layers such 
that the interlayer coupling energy in fact vanishes at the 
mean field level. 

For Cs2CuBr 2 (J/ J' = 2 and D — 0) series expansions 
predict M ss 0.17. 17 It would be interesting to make 
precise neutron scattering measurements of the magni- 
tude of the ordered spin moment in both CS2CUCI4 and 
Cs2CuBr4 to compare directly with the theoretical cal- 
culations presented here. 



V. MAGNON DISPERSION 

A. Cs 2 CuCl 4 

In this subsection we present series expansion results 
for the magnon dispersion for parameters relevant to 
CS2CUCI4, discuss the effects of the DM interaction on 
this dispersion, and compare it to the dispersion obtained 
from spin- wave theory with 1/S corrections 9 ' — and to the 
experimental dispersion obtained from inelastic neutron 
scattering with unpolarized neutrons, 3 



TABLE I: Comparison of values for the incommensuration e 
and sublattice magnetization M for CS2CUCI4 obtained from 
experiments and various theoretical approaches. Unless noted 
otherwise, the values of the exchange interactions J and J' 
used in the theoretical calculations are J — 0.374 meV and 
J' = 0.128 meV, corresponding to a ratio J /J' — 2.92. D 
denotes the strength of the Dzyaloshinskii-Moriya interaction. 
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FIG. 3: (Color online) Reciprocal space diagram for the 
triangular lattice. The path ABCOAD (A = (it,0), B = 
(tv,ty/V3), C = (27r/3,27r/V3), O = (0,0), and D = (2w,0)) 
and the three sets of points denoted (l)-(3) are cuts along 
which magnon dispersions are plotted in subsequent figures 
in the paper. 



In Fig. H we show series expansion results for the 
magnon dispersion along the fc-space path ABCOAD in 
Fig. [3] The dashed blue curve is for J / J' = 2.92 and 
D = 0, while the full red curve includes the effect of a 
finite D/J' = 0.16. We note the following features: 

(i) Along AB, which is perpendicular to the chain di- 
rection (6 axis), the excitation energy for D/J' = 0.16 
is significantly enhanced (by more than a factor of two) 
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FIG. 4: (Color online) Magnon dispersion uik, as calculated 
from series expansions, for parameters relevant to CS2CUCI4, 
plotted along the path ABCOAD shown in Fig. |3] When the 
DM interaction is included the magnon energy is significantly 
enhanced along AB and near the point D, and a gap opens 
up at the ordering wavevector k = Q. 



over that for D = 0. The dispersion along AB is very flat 
in both cases. 

(ii) The D 7^ dispersion has a gap at k = Q while the 
D = dispersion is gapless there. This follows from sym- 
metry arguments. For D = the model has full SU(2) 
spin rotation symmetry. The helical order then leads to 
gapless excitations (Goldstone modes) at k — and at 
the ordering vector k = Q . The Goldstone mode at 
k — Q is associated with rotations of the ordering plane 
of the spins, which is arbitrary when D = 0. The SU(2) 
symmetry is broken by the the DM interaction which 
fixes the ordering plane to coincide with the plane of the 
triangular lattice. Thus rotations of the ordering plane 
costs a finite energy for D ^ 0, which creates a gap in the 
magnon dispersion at k = Q^ i 25 i 26 i 27 (it is important to 
note, however, that to get gapless excitations at the ap- 
propriate k vectors in the series calculations we need to 
bias the analysis as discussed in some detail in Ref. [23|; 
an unbiased analysis always gives a gap.) 

(hi) Near the point D the excitation energy for D/J' — 
0.16 is significantly enhanced (by approximately a factor 
of three) over the D = case. 

(iv) Overall, the high-energy parts of the dispersion are 
quite insensitive to the DM interaction, while the low- 
energy parts are quite sensitive. A notable exception to 
the latter is the region around fc = 0, since the Goldstone 
theorem dictates gapless excitations at k = regardless 
of whether D is zero or not (this is because the k = 
Goldstone mode is associated with long- wavelength rota- 
tions of the spins within their ordering plan&2 i 25 ' 26 i 27 ). 
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FIG. 5: (Color online) Magnon dispersion for J/J' = 2.92 and 
D/J = 0.16 (solid points from series, black solid line from 
spin- wave theory with 1/S correction) compared to experi- 
mental dispersion in CS2CUCI4 (squares from Ref. 0, dashed 
line is experimental parameterization) along the path AB- 
COAD in Fig. [3] Compared to the dispersion from spin- 
wave theory, the series dispersion is enhanced along the J 
bonds and decreased perpendicular to them. Note that theo- 
retical and experimental dispersions cannot be directly com- 
pared near the ordering wavevector k — Q, as the gap in 
expected at this wavevector cannot be resolved in the unpo- 
larized neutron scattering experiments (see text about how 
the experimental dispersion is extracted). 



2. Comparison with spin-wave theory 

Next we compare our theoretical dispersion for D/J' = 
0.16 with that obtained from spin-wave theory with 1/S* 
corrections (LSWT + 1/5 for short) £±2 Both disper- 
sions are plotted in Fig. [5] (black dots and full line, 
respectively). We see that compared to the spin- wave 
prediction, the excitation energy is increased in the b di- 
rection along which neighboring spins are coupled by the 
strong J bonds, and decreased in the perpendicular c 
direction. This corresponds to an upwards renormaliza- 
tion of J and downwards renormalization of J' with re- 
spect to the spin-wave prediction, which thus effectively 
makes the system appear more one-dimensional. (The 
renormalizations of J and J' with respect to linear spin- 
wave theory (LSWT) are even bigger^ as the LSWT + 
1/5 dispersion itself is in turn characterized by an up- 
wards (downwards) renormalization of J (J') compared 
to LSWT.) The dependence of the magnon energy on 
the c-component of the wavevector is weak, giving the 
dispersion a pronounced one-dimensional character. 



3. Comparison with experimental dispersion 

In Fig. [5] we also show experimental results for 
the magnon energies obtained from inelastic neutron 



scattering. 3 The open symbols are data points corre- 
sponding to the positions of the strongest peaks in the 
unpolarized neutron scattering data. The red dashed line 
is a fit of these points to a dispersion given by 



<-Ok = \ (Jk - Jq) (JQ+k + JQ-fe)/2 - Jq 



(12) 



where Jk = J cos kb + 2J'cos(fc b /2)cos(v / 3fcc/2). The 
functional form of Eq. (fT2|) is identical to the expres- 
sion from linear spin-wave (LSWT) theory in the ab- 
sence of the DM interaction, but in Eq. (|12[) the bare 
exchange parameters J and J' that would be used in 
LSWT have been replaced by effective, renormalized pa- 
rameters (J and J', respectively) whose values are cho- 
sen to give the best fit to the experimental data points. 
Hence the ratios of these renormalized values to the bare 
values give a measure of the strength of the quantum 
fluctuations in the system. In Fig. [5] J — 0.61(1) meV 
and J 1 = 0.107(10) meV£ 

Before comparing the theoretical and experimental 
curves, however, we first discuss the interpretation of the 
experimental data in some more detail, i.e. why for a 
given k the location in energy of the strongest peak in 
the data can in most cases be identified with the magnon 
energy Uk- To this end, we consider the differential cross 
section for inelastic scattering of unpolarized neutrons, 
which is proportional to2£ 



sm 



9 k S aa {k, w) + (1 + cos" 9 k ) S bb (k, w) 



(13) 



where 9k is the angle between the scattering wave vec- 
tor k and the axis (a) perpendicular to the ordering (be) 
plane of the spins, and S aa (k, to) and S bb (k, to) are diago- 
nal components of the dynamical structure factor defined 
as (/i, v = a,b, c) 



1 r°° 

S^(k,to) = — dt V(5 /i (0,0)5 !y (i?,t))e- l(fc - R - 
2ir I z — ' 

J -co R 



A) 



(14) 



When k is in the be plane (as in our calculations), 9k 
7r/2, in which case Eq. (fT3|) reduces to 



S aa (k,u;) + S bb (k,Lu) 



(15) 



Spin- wave theory predicts sharp one-magnon peaks in the 
out-of-plane correlations S aa (k, u>) at ui = u>k (referred to 
as the principal mode) and in the in-plane correlations 
S bb (k, ui) at both u> — LOk+Q = ^t, anc ^ 10 — U k-Q = ^k 
(referred to as the two secondary modes), where u>k is the 
magnon dispersion. Hence if the spin- wave prediction for 
Eq. (fl~5|) were plotted as a function of to for fixed k, three 
peaks would be observed, at to — u>k and wj^. Although 
spin-wave theory does not give the correct fc-dependence 
of the three modes, it is expected that the three-peak 
structure it predicts (for a system with helical magnetic 
order) is qualitatively correct. In Fig. [S]we have plotted 
this structure for k varying along the path in Fig. [3] 
(the energy of the modes was calculated from spin-wave 
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theory with 1/5 corrections^) . For example, at point B 
the spin-wave calculation predicts that lo^ lies well below 
the two other modes in energy. This principal mode is 
also predicted to have the strongest intensity at B. Thus 
it appears quite safe to assume that the lowest-energy, 
strongest-intensity peak in the experimental spectrum at 
that fc-vector should be assigned to uj k . At high energies 
the intensity of the out-of-plane correlations S aa is also 
predicted to be bigger which justifies an assignment of the 
location of the strongest peak to u>k at high energies. On 
the other hand, at low energies close to k = Q, where the 
primary mode u>k is expected to be gapped and located 
inside the V-shape of the wjT mode (see Fig. [6]) , the S aa 
part can not be separately resolved. The experimental 
data indicate that there is no (or very little) overall gap 
at k — Q. This is what one would expect, because for 
a spiral with finite DM interaction one of the in-plane 
modes (u>^) would still be gapless. 




(7t,0) (n,m ) (271/3,271/3") 



(0.0) 



(71,0) 



(2ti,0) 



FIG. 6: (Color online) The magnon modes cuk (principal 
mode) and uik±Q = wj^ (secondary modes), as calculated 
from spin- wave theory with 1/S corrections^ for CS2CUCI4 
(J/ J' = 2.92 and D/J' = 0.16), plotted along the path AB- 
COAD shown in Fig. [3] 



For point B, uj^ is quite unambiguously determined, as 
discussed above, and this point is sensitive to the pres- 
ence of the DM interaction; the theoretically calculated 
energy increases with D. This makes the interchain dis- 
persion in better agreement with the experimental data 
when D is included, although there is still a slight over- 
estimate of ~ 20%. Along the AD direction the agree- 
ment between experiments and series is essentially per- 
fect, while along the chain direction OA, where the dis- 
persion is maximal, the theory makes a slight underes- 
timate of ~ 10%. The experimental dispersion relation 
was also measured along off-symmetry directions in the 
2D Brillouin zone and Fig. [7] shows that the series results 
(black circles) compare very well with experimental data 
(red squares) for those directions too. 
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FIG. 7: (Color online) Comparison of the experimentally mea- 
sured dispersion relation in CS2CUCI4 (red squares) and se- 
ries expansion results for J = 0.385 meV, J / J = 2.92 and 
D/J' = 0.16 (black circles), along paths (1), (2) and (3) in 
Fig. [3] Calculations were made at the same (k x ,k y ) where 
experimental dispersion points were measured. 




FIG. 8: (Color online) Magnon dispersion for Cs2CuBr4 (J = 
2 J' and D — 0), calculated from series expansions. Compared 
to the dispersion from linear spin-wave theory (red dashed 
line, Eq. (JT2J ) the series dispersion is enhanced along the J 
bonds and decreased perpendicular to them. The dispersion 
looks qualitatively rather similar to the case J/ J' = 2.92, 
D = in Fig. [4] although in that case the energy is practi- 
cally flat along AB (i.e., perpendicular to the chains) while 
here the energy varies somewhat along AB, giving a less one- 
dimensional dispersion, as expected. The horizontal black 
dashed line is at uo/ J — it/2 which is the maximum energy in 
the one- dimensional limit J/ J' =00. 
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B. Cs 2 CuBr 4 

In Fig. [8] we show series expansion results (black sym- 
bols) for the magnon dispersion for J/ J' = 2, the ratio 
appropriate for Cs2CuBr 4 , along the fc-space path AB- 
COAD in Fig. [3] The dispersion obtained from linear 
spin-wave theory, Eq. (|12[) is also shown for compari- 
son. We note the following features: (i) With respect 
to LSWT the series dispersion is renormalized upwards 
in the direction parallel to the chains (i.e., in the di- 
rection along which the exchange constant is strongest) 
and downwards in the perpendicular direction. A simi- 
lar kind of upward renormalization occurs both for the 
square lattice (J = 0) and decoupled chains (J' = 0). (ii) 
Along the AB direction perpendicular to the chains, the 
series dispersion has flattened and the excitation energy 
has decreased significantly compared to the case J = J' 
case (Fig. 1 in Ref. |22( ). Features (i) and (ii) persist as 
J I J' is increased further to 3, as seen in the D — curve 
in Fig. QJ which except for the more pronounced flatness 
along AB looks qualitatively very similar to the J / J' = 2 
dispersion in Fig. [H 

VI. CONCLUSIONS 

In this paper we have presented series expansion calcu- 
lations for various ground state properties (ground state 
energy, ordering wavevector, and sublattice magnetiza- 
tion), and for the magnon dispersion, in the helically or- 
dered state of spin-1/2 frustrated antiferromagnets on an 
anisotropic triangular lattice. For the parameter values 
considered here the model Hamiltonians we have studied 
are expected to give a good description of the magnetic 
properties of Cs2CuBr4 and CS2CUCI4. 

For parameters appropriate for Cs2CuBr4 (J/ J' = 2) 
we have made specific predictions for large and non- 
uniform renormalizations of the dispersion relation com- 
pared to the classical (linear spin-wave theory) result (see 
Fig. [5]). There is an upward renormalization along the 
chains and a downward renormalization along the inter- 
chain direction. These predictions for Cs2CuBr4 should 
be testable by future neutron scattering experiments. 

For CS2CUCI4 we have compared our theoretical pre- 
dictions for the dispersion with experimental results 
from neutron scattering experiments and with spin- wave 
theory that includes 1/S corrections. The calculated 
magnon dispersion shows large quantum renormaliza- 
tions compared to the classical result; as for Cs2CuBr4 
the renormalization is upward along the chain direc- 
tion and downward along the interchain direction. Thus 
quantum fluctuations make the dispersion look more one- 
dimensional. These predictions from series are in good 
quantitative agreement with the experimental observa- 
tions. In contrast, spin- wave theory with 1/5 corrections 



predicts considerably smaller renormalizations and thus 
underestimates the effects of quantum fluctuations on the 
dispersion. 

The agreement between theory and experiment im- 
proves further when the Dzyaloshinskii-Moriya (DM) in- 
teraction is taken into account. The DM interaction is 
found to significantly increase the incommensuration of 
the ordering wavevector with respect to the antiferro- 
magnetic wavevector in the chain (b) direction. Another 
important consequence of the DM interaction is that it 
opens up a gap in the magnon dispersion at the order- 
ing wavevector. We have made specific predictions for 
the size of this energy gap at the ordering wavevector 
to longitudinally-polarized excitations (along the direc- 
tion normal to the spiral plane) which could be tested by 
future high-resolution polarized neutron scattering exper- 
iments. 

Another very interesting issue, at least from a theoret- 
ical point of view, is the role of the DM interaction in 
CS2CUCI4 in stabilizing a magnetically ordered ground 
state. We find here that inclusion of the DM interaction 
strengthens the magnetic ordering. In fact, several re- 
cent paper a 10 ! 11 ' 12 ! 13 have predicted that in the absence 
of the DM interaction the ground state is not magneti- 
cally ordered. The proposed alternatives to a magneti- 
cally ordered phase for D = include a spin liquid ("al- 
gebraic vortex liquid")^ and a state with "staggered" 
dimerization^ The latter state could be studied using 
series expansions; it would be interesting to compare its 
energy with the helically ordered phase. The fact that a 
different dimerized state (with "diagonal" dimerization; 
see Fig. 4a in Ref. flTt ) that we have studied with series 
expansions is extremely close in energy to the helically or- 
dered state (for a large regime of J/ J' including J/J' « 3 
we find that the energies of the two phases lie virtually 
on top of each other) is strongly suggestive that the state 
with staggered dimerization proposed in Ref. [l3| may be 
highly competitive as well. 
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TABLE II: Series coefficients for the magnon dispersion for the isotropic triangular-lattice model, calculated for J = 2.92, 
J' = 1, q = 1.64, D = 0.16, and t = 4 in Eqs. @ and (|10[1 . Nonzero coefficients c rimi „ in Eq. pip up to order r — 8 are listed. 
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